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Abstract 

The (M, _fC)-reduced non-autonomous discrete KP equation is linearised on the 
Picard group of an algebraic curve. As an application, we construct theta func- 
tion solutions to the initial value problem of some special discrete KP equation. 

1 Introduction 

The non-autonomous discrete KP equation (ndKP) is given by the formula [6]: 

(b(m) - c(n)) • tit}jL+i.n+i + (c(n) - a(t)) • /^ +1>n /^+i 
+ (a(t) - b(m)) ■ f m , l+l f+ + \ n = 0, t, m, n e Z. 

With constraints a(t) = 0, b(m) = 1, c(n) = 1 + 8 n and „ = /^J*R- n , the 
ndKP reduces into the following form: (f n K ~ mM ; = /m n); 



J n+1 /n+1 /-. | r \ J nJ n+1 c Jn /n+1 

rt+M+K ft \ + n+1 > ft ft+K ~ ~°n+l t+K ft+M +h, 

/n+1 /n+1 /n+l/™ /« /n+1 



j-t jt+K jt+K jt-M+K 

Define 7* := (1 + 5 n+1 ) • " and 7 n * := 8 n+1 ■ n t+K n t ^ M+K ■ Then we 

/n+l/n _ /n+1 /« 

derive the following discrete system: Vn, i G Z, 

i n — i n-l + v n "n-l> ^-v 

jt-Myt-K 

-*n 

The system (|l.llll.2|) is called (M, K) -reduced non- autonomous discrete KP 
equation (rndKP). The term 'non-autonomous' derives from the freedom in 
the parameters S n . If we assume an extra constraint 6± = 6% = 63 = • • • , this 
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system reduces to an autonomous system ((M, K) -reduced autonomous discrete 
KP equation (rdKP)). 

In this article, we study the rndKP with the periodic boundary condition: 



■N' 



V 

1 n 



yt 

V n+N> 



N e N. 



(1.3) 



The present paper is a generalisation of the method to solve the generalised 
periodic discrete Toda equation introduced in the papers [TJ [2] . We show here 
that this method is also applicable to the quite general case of the rndKP and 
prove a linearisation theorem (theorem 12 . 1 2[) , which illustrates the geometric 
information of the discrete system. 

In some special situation, theta function solutions of the initial value problem 
are constructed. In section [31 we derive an explicit formula for the solutions of 
the rdKP, which is a reduction of the rndKP. 

Important remark We can assume g.c.d.(M, K) = 1 without loss of gen- 
erality. (See (fTTTjl . (fTT2)) .) Aside from this, we assume g.c.d.(M +K,N) = 1 
in Sections [H and [3] by technical reason. The general cases will be discussed in 
Section [U 

Notation: For a meromorphic function / over a complete curve C, (/)o 
(resp. (/)oo) denotes the divisor of zeros (resp. poles) of /. Let (/) := (/) — 
(/)oo- Div d (C) means the set of divisors over C of degree d and Pic d (C) 
means the quotient set defined by Pic d (C) = Div d (C) /(linearly equivalent). 
For an element V £ Div d (C), [D] means the image of T> under the natural map 
Div d (C) ^Pic d (C). 



2 Inverse scattering method 

The rndKP equation (|l.lH1.3p has the following matrix form: 

L t{y)Rt{y) = R t -Ai(y)Lt- K (y), 

where 

(vi i \ /it i 



\ y 



i 



Rt(y) 



\y 



(2-1) 



i 



and y is a complex parameter. Let 

X t (y) ■= L t _ {K _ 1)M {y) ■ ■ ■ L t -2M{y)Lt-M(y)L t (y)x 

R t (y)R t -K(y)Rt-2K(y) ■ ■■Rt-{M-i)K{y), 

then (|2.ip becomes 

X t (y)Rt-Mi<(y) = R t -MK(y)X t -K(y), 



(2.2) 
(2.3) 
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or equivalently 

L t -MK(y)X t (y) = X t - M (y)Lt-MK(y)- (2.4) 

Because M and K are co-prime, the characteristic polynomial of X t (y) does not 
depend on t. Let C := {(a;, y) G C 2 | det - i£) = 0}. Of course, C is 

also independent from i. We call the completion C of C the spectral curve of 
the rndKP equation. 

Remark 2.1 5j/ applying l\2.1\i repeatedly, we can transform (|£.£?p mio 

= Rt-MK {y)Rt-(M+i)K(y) ■ ■ ■ Rt-(2M-i)K(v)* 

L t -(2N-l)M(y) ■ ■ ■ L t -(K+l)M(y)L t -KM(y)- (2.5) 

2.1 Properties of the spectral curve 

As a starter, we list some fundamental properties of the spectral curve C in this 
section. In the rest of this article, we always assume C to be smooth. Moreover, 
we also assume that g.c.d.(M + K, N) — 1 in Sections and [3] unless otherwise 
is stated. 

Denote the set of N x N matrices by Mjy(C) and the subset of diagonal 
matrices by V C M N (C). For a matrix X G M N (C) and subsets A, B C M N (C), 
let A+X := {a+X | a G A}, AX := {aX \aE A}, A+B := {a+b \ae A,be B} 
and AB := {afr | a G A, 6 G B}. 

f° 1 \ 



Let S be the N x N matrix 5 







'•• 1 
V?/ 0/ 

The polynomial detpQ(y) — xE) is of degree AT w.r.t. x, and of degree 
M + K w.r.t. y. Then the projection p x : C 3 (x, y) h-> x G P is (M + K) : 1, 
and the projection p„ : C 9 (a;, y) i— ► y G P is iV : 1. 

Let U S := (nf=i^ +fc ) ' (Il^i^ +k ). (J'6 {L2,...,iV}). By the 
quantity C/j is invariant under the time evolution. 

Proposition 2.1 If g.c.d.(M +K, N) = 1, the curve C has the following special 
points: 

(i) M points Aj : (x, y) = (o, ULi I n 3K ) , J = 0, 1, . . . , Af - 1 . 

(ii) K points Bj : (x, y) = (o, nLi Vr7 jM ) , 3 = 0, 1, • • • , K - 1. 



(hi) iV poinis Q 3 : (x, y) = ([/,-, 0), j = 1,2, . . . , N. 
(iv) a unique point P : (x, y) = (oo, oo). 
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Proof. Let (0, y) G C. Then we easily derive 

K-l M— 1 

H dct {Lt-iM) • n det (R t - jK (y)) = 0, 

3=0 j=0 

which implies (i) and (ii). Part (iii) follows from the fact that L t (0) and i?t(0) 
are upper triangular. 

(iv): For a point (x,y) G C, there exists a non-zero TV-vector v(x, y) such 
that X t (y)v(x,y) = x ■ v(x,y). Because the matrix X t (y) is contained in the 
subset (r + S) M+K = T + TS + ■ ■ ■ + TS M+K - 1 + S M+K , it follows that 

(7o + 7i 5 + ' ' ' + 7m+a'-i5 m+k - 1 + S M+K ) -v = x-v, (2.6) 

where 70, 71, ... , Jm+k-i are diagonal matrices. 

Define a new parameter k by y — k~ N which is assumed to be zero near 
P. Let £/v be a N-th. primitive root of unity. For all j € {0,1, . . . ,N — 1}, the 
vector 

v :=(((ik) N -\(<;ik) N - 2 ,...,(<;ik), if 

satisfies the formula: S ■ Vq = (Cjv^) 1 ' v o- Then, from (|2.6p we obtain 
(Cat^) - ' ^0 = x ■ Vo + (higher term) near k = 0, which implies (x,y) ~ 
((( J N k)- M ~ K , k~ N ) when (ar, y) G C tends to infinity. Because M + if and N 
are relatively prime, we can choose an appropriate branch of k around a unique 
point P such that 

X = k^W + • • • , y = k - N + ■■■ , 

near P. ■ 

From the proof of proposition 12.11 one obtains more detailed information on 
the point P G C. 

Corollary 2.2 There exists a local coordinate k around P such that 
x = fc-< M + A ") + . . . , y = k~ N + ■■■ . 



Corollary 2.3 Let X t (y)v(x, y) — x ■ v(x, y). Then, around P , it follows that 
v(x, y) ~ (k 1 ^^ 1 , k N ~ 2 , . . . , k, 1) T (up to a constant multiple). 



Proof. Because g.c.d.(M + K, N) = 1, the solution of the vector equation 
is expressed as v(x, y) = (k 1 ^^ 1 , k N ~ 2 , . . . , k, 1) T + (higher) up to a constant 
multiple. ■ 
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Remark 2.2 The proof of vrovosition \2. 1\ (i). (ii) implies that the set 

zs invariant under the time evolution. It then follows that {Il^n'II^ra} = 
{n^ +K >n-^ +M }- To avoid a non-interesting solution J*+ M = V^, F„ t+Ar = 
^n+i °/ ^ e rn dKP \1.1W1.^ , we should assume the extra constraint Yi n I n +M — 
Yin 7^ Yin Vn +K ~ Yin ^ri ^ n addition to the rndKP. In fact, this constraint 
is enough to guarantee the existence of the unique solution. See section \2.3[ 

Next we consider the behaviour of Qj (J = 1,2, ... , N). The position of Qj 
is invariant under the time evolution. In this paper, we restrict ourselves to the 
following two typical cases: 

(a) All Qj are distinct. (b) Q\ = Q2 = ■ • • = Qn- 

Note that in the case (b), the system (jl.HHl.3P reduces to the rdKP. 

In the case (a) 

The equation X t {y)v(x, y) = x ■ v(x, y) becomes 

Z t (Q) • v = Uj ■ v, at Qj. 
Because X t (0) is upper triangular, the eigenvector v takes the form 

v = {d l7 d 2 ,...,d ]7 0,...,0) T , dj^O. (2.7) 

In the case (b) 

Let Q := Q\(= Qi = ■ ■ ■ — Qn)- Arguments similar to those in the proofs of 
corollarv l2.21 12.31 prove the following: 

Proposition 2.4 If all the points Qj coincide, there exists a local coordinate k 
around Q = Q\ such that 

x = [/i+--- , y = k N + ■ ■ ■ , 

and the eigenvector X t (y)v(x,y) = x ■ v(x,y) satisfies 

v(x, y) ~ (1, k, k 2 , . . . , k N ~ 1 ) T , (up to a constant multiple). 

m 

Let v(x,y) = (gi{x,y), . . . , gN(x,y)) T be an A-vector function (defined up 
to a constant multiples) such that X t (y) ■ v(x,y) — x ■ v(x,y). By the above 
arguments, we have: 
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Proposition 2.5 The meromorphic function gj/gj + i, (j < j + I < N) has: 

(i) I zeros at P 

(ii) at least one pole at Qj . 

Define the divisors T>\ and T> 2 to be minimal positive divisors on C such 
that 

{93/ 9n) +Vi> -(Qj + Qj+i + ■■■ + Qn-i), Vj, (2.8) 
( gj /9i)+V 2 >-(j-l)P, Vj. (2.9) 

These divisors were first studied in [3], where it has proved that T>\, T> 2 are 
general and degVi — degT> 2 — genus(C). 

2.2 the eigenvector mapping 

Let p be a point on a smooth curve C and k be a local coordinate around 
p. For a meromorphic function /, ord/(p) denotes the largest integer r such 
that lim^p \k~ r f(q)\ < +00. For a vector function v(p) = (fi)i, we define 
ordv(p) := min, [ord (/i(p))]. 

An isolevel set Tc is the set of matrices X(y) (eq. (|2.2p ) associated with the 
spectral curve C. Let <? := genus(C). Now we construct a map from Tc to 
Pic 9+Ar_1 (C) called the eigenvector mapping. 

Let X = X(y) be an element of Tc- If (^,y) G C, there exists a complex 
A^-vector v(x,y) such that X(y)v(x,?/) = xv(x,y), up to a constant multiple. 
Then there exists a Zariski open subset C° of C over which the morphism 
C° 9 (a;,y) 1— * v(x,y) 6 p^- 1 j s uniquely determined. Moreover, for a smooth 
C, this morphism can be extended uniquely over the whole of C. Denote this 
morphism by ^> x ■ C — > P^ -1 . 

The eigenvector mapping (pc : Tc Pic d (C) (d = g + N — 1) is a map 
defined by the formula: 

O c (<Pc(X)) = V'xiOrN-xil)), (2.10) 

where Opjv-i(l) is the invertible sheaf of hypcrplane sections over P^ -1 . Note 
that it is nontrivial to prove that ipc(X) € Pic d (C) (see [I] §2). 

Let (Xi : X 2 : • • • : Xn) be the homogeneous coordinate of P^ -1 . The eigen- 
vector mapping illustrates the geometric interpretation of the general divisors 
T>i and T>2 ( section l2.ip . In fact, (|2.8p implies that 

T> 1 + Q 1 + Q 2 -\ h Qjv-i is the pull-back of {X N = 0}, 

and (|2.9p says 

V 2 + (N - 1) • P is the pull-back of {X 1 = 0}. 

These facts imply the following: ipc(X(y)) = [T>i + Qi + Q 2 + ■ ■ ■ + Qjv-i] = 
[V 2 + (N-1)-P]. 

Let l)(X(y)) := T> 2 . This divisor will play an important role for constructing 
a tau function solution of rdKP. See the next section. 
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Remark 2.3 Because T> x + Q 1 + Q 2 H h Qn-i and T> 2 + (N - 1) • P are 

linearly equivalent to each other we have 

{91/ 9n) = T> (X(y)) + (N - 1) • P - X>! - (Q 1 + Q 2 + ■ ■ ■ + Q N ^). (2.11) 

Remark 2.4 Let X(y) ■ v(p) = x ■ v(p) and p — (x, y) <E C . Equation \2. 10$ is 
equivalent to ipc(X) = —^2 peC (ordv(p))-p . 



The following theorem is essentially obtained in van Moerbeke, Mumford [I] . 
Theorem 2.6 The eigenvector mapping ipc ■ Tc ~ * Pic d (C) is an embedding. 

2.3 shift operators 

Consider the N x N matrix Xt(y) defined by 1|2.2|) and the associated spectral 
curve C. Let <r, \ik and (j,m be the isomorphisms on Tc defined by: 

a(X t (y)) := SX t (y)S-\ (2.12) 
li K {X t {y)) := R t _ {M _ 1)K (y) ■ X t (y) ■ {R t -(M-i) K (y)}-\ (2.13) 
H- M (X t (y)) := L t _ MK (y) ■ X t (y) ■ {L^ M k{v)}~\ (2.14) 

where S is the matrix defined in section [2~T1 By (|2.3H2.4p . we have fiK(Xt) = 
X t+ K and ii-M{X t ) = X t -u- For the rndKP (|1.1H1.3[) . a is the n-shift op- 
erator: n 1 y n + 1 and /j,k and (A—m are the t-shift operators: t 1— > t + K, 
1 1 — ^ t — M. Because X and M are co-prime, an appropriate combination of \xk 
and /1 Af defines the unit time evolution 1 1— » t + 1. 

We start with the linear problem: 

X t (y) -v(x,y) = x-v(x,y), v(x,y) = {gi{x,y))f =1 . (2.15) 

This linear equation is decomposed into the following infinite dimensional form: 
for an infinite vector (<7i)i g z, 

a i.O ' 9i + a i,l ' 9i+l + ' ' ' + di.M+K ' 9i+M+K = x ' 9i, ( a i+N.j = a i.j) 

(2.16) 

gi+N = y-9i- (2.17) 

The matrix equation (|2. 15[) can be interpreted " (|2.16[) with constraint (|2.17p ". 
However, interchanging the roles of these two equations, i.e., interpreting " (I2.17P 
with constraint (|2.16p " , we arrive at another matrix equation: 

Y t {x) w =y-w, where w = { gi )fi+ K . (2.18) 
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ai a 2 1 \ / gi \ / gi 

Example 2.7 For an equation | y b\ hi I \ gi I = x \ gi ) , the asso- 

c 2 y y ex) \g 3 J \g 3i 

ciated new matrix equation is: 



b 2 (ai - x) 



a 2 b 2 - bi + x 



(ai - x)(ci - x - b 2 c 2 ) a 2 (ci — x) — c 2 (a 2 b 2 — bx — x) 



= y 



We call the linear problem (|2.15|) the x-form and the linear problem (|2.1£ 
the y-form. 



2.3.1 shift operators and the x-form 



Due to (|2.12| - |2~h3| . the shift operators a and fi act on the eigenvector of the 
i-form equation (|2. 15[) by: 

a : v S v, fi K ■ v ^ {Rt-(M-i)i<} v, (J,-m ■ v ^ {Li-mk} v. 
The following lemma is easily proved: 

Lemma 2.8 det S = { 7 }) N+1 y, det R t - {M -i)K = Tin J*- (M - 1)K - - 
det L t - 



t-MK 



n v 
1 in ' 1 



t-MK 



2.3.2 shift operators and the y-form 

The y-form representation of the shift operators a, [Ik, fJ>—M are more com- 
plicated. Let Ei := -(01,0 - x)/cli i m+k, E 2 := —ai,i/ai,M+K, ■ ■ ■ ,E m +k ■= 
-ai,M+K-i/ai,M+K- Then (|2.16|) becomes g M +K+i = E i9i- 
Define three new matrices S* , R* and L* by: 



S* 



R* 



( 



1 



E 2 



\Ei E 2 



L* 



\ 



1 

Em+k-i Em+k J 



± M+K-l 

Em+k-i 



I M+K + Em+K J 
\ 



\ E\ E 2 



*M+K-X 

Em+k -i 



M+K 



+ Em+k / 



(2.19) 



(2.20) 



(2.21) 
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where I~ = f-{M-l)K &nd y _ = yt-MK The matrices 5 * ; R * and j* are 
the y-form version of the matrices S, Rt-(M-i)K and L t -MK "with constraint 
9m+k+i = Eigi" . The shift operators a and \i act on the eigenvector of the 
y-form equation (|2.18|) by: 

(7 : w i — ► S*w, fix ■ w i ► _R*k;, M-a/ : w l_ * L*n>. 

Lemma 2.9 detS* = (-l) M + K ■ (U x - x), detiT = detL* = {-l) M+K+1 -x. 

Proof. The calculation is cumbersome but elementary. We will prove this 
lemma in Appendix. ■ 



2.3.3 geometric interpretation of x-form and y-form 

Consider the projections p x : C — > P and p y : C — * P (section l2.1j) . Recall that 
p x is (M + K) : 1 and p y is N : 1. Denote J 7 := O c (tpc(X t )). 

Because the a;-form representations of cr, are independent from x (section 
I2.3.ip . for fixed y e P and its pre-image Py X {y) = {(x\,y), . . . ,(xN,y)}, the 
matrices S, Rt-jM-i)K and Lt-MK act on the vectors v(xi,y), . . . ,v(xjsr,y) 
simultaneously. Q| □ 

On the other hand, for generic y, the vectors v{x\, y), . . . , v(xn, y) should be 
linearly independent because they are eigenvectors belonging to distinct eigen- 
values. What happens if we choose y such that det S(y) = ? This seemingly 
leads to a contradiction, if one believes the the column vectors of the singu- 
lar matrix S(y) ■ (v(x±, y), . . . , v(xn, y)) are linearly independent. However, 
realizing the fact that the eigenvectors are only determined up to a constant, 
this problem is easily solved and we conclude that J2i=i or d (S(y)v(xi, y)) > 
Eiliord (v(x h y)). 

More precisely, the statement of lemma 12.81 can be interpreted as follows: 

• J2iLi 0ld ( Sv ( x i^)) = Y^iLi old ( v ( x i^)) + 1 > where (x t ,0) G C, i = 
1,2,. ..,N. 

. Lcty :=n„^r (M_1) ^- Then 

Eili ord {Rv(xi,y )) = Eili ord ( v ( x i> yo)) + L (xi,yo) e C. 

• Let yi :=Un V n~ MK - Then 

J2fLiOTd(Lv(x t , yi )) = E^IiOrd(v(a; l ,yi)) + 1, (x t ,yi) G C. 

If X(y) has an eigenvalue x' of multiplicity m > 1, we should choose the vectors 
v(p), v'(p), . . . , -i/ m_1 ) (p), where (p) is the fc-th differential of v with respect to the local 
coordinate around p = (x',y). 

2 Geometrically, this means that S, R and L act on the push-forward {p y ) t T. 
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Similar arguments in the case of the y-form representations yield the follow- 
ing form of lemma 12.91 

• EZt 1 ^(S*w(U 1 ,y i ))=Y,Zt 1 ord(w(U 1 ,y i )) + l, (E/i, W ) e C. 

• E -It' ord (R*w(0, jft)) = Efit 1 OTd (*"(°> W)) + 1, (0, Vi) 6 C. 

• E^f ord (L*t»(0, = ES 1 ord (™(0, + 1, (0, y 4 ) G C. 

Combining these data, we obtain the following proposition: 
Proposition 2.10 Let Q\ : (x,y) — (U\,0) and 

A 3 :(x, y) = (0, n„ t~ {M ' 1)K ), -jK = t - (M - l)K (mod M), 
Bi :(x, y) = (0, U n V*~ MK ), -%M = t— MK (mod K) 

(proposition \2.1\ . Then. 

(i) ord {Sv(Qi)) = ord (»(Q X )) + 1, (ii) ord (Jto^)) = ord (u^)) + 1, (iii) 
ord (Lv(Bi)) = ord (« (B,)) + 1. 

Proof. We prove (i). By construction of the x-form and the y-form, we have 

ord (Sv(p)) = ord (v(p)) + 1 ord (S*w(p)) = ord (w(p)) + 1. 

On the other hand, because a regular matrix is invertible, 

det S(y) 7^ 0, oo or det S*(x) ^ 0, oo ord (Sv(x, y)) = ord (i>(x, y)). 

These facts prove proposition (i). Clearly, similar arguments will prove (ii) and 
(iii). ■ 

2.3.4 shift operator at the infinity point 

At P, the actions v(P) h-> Sv(P), v(P) ^ Rv(P) and v(P) ^ Lv(P) are 
directly computable. 

Proposition 2.11 (i) ord (5u(P)) = ord(u(P)) - 1, 

(ii) ord (Rv(P)) = ord (v(P)) - 1, (iii) ord [Lv{P)) = ord (v(P)) - 1. 

Proof. The Proposition is readily proved by Corollaries 12.21 and 12.31 ■ 
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2.4 linearisation theorem 

From the above calculations, we obtain the linearisation theorem representing 
the flow of the rndKP equation on the Picard group of the spectral curve. 

Theorem 2.12 (I): Let T> be the divisor T> = P — Q\, Then the following 
diagram is commutative. 

T c -> Pic d {C) 
a I I +p] . 

T c -> Pic d {C) 

(II) x ; LetSj (j = 0, 1, . . . ,M-1) fee tfie divisor £j =P-Aj andt= -(j + l)K 
(mod M). The following diagram is commutative. 

T c - Pic d {C) 

M-K" 4 I +&] ■ 

r c - Pic d {c) 

(II) 2 : Lei J^j (j = 0, 1, . . . , K - 1) be the divisor Fj = P - Bj, and t = -jM 
(mod K) . The following diagram is commutative. 

T c Pic d {C) 

M-m I I + ■ 

T c -» Rc d (C) 

Proof. The theorem follows immediately from Remark 12.41 and Proposition 
l2~T01 ■ 

We should note the fact that the position of the points Qj (proposition 1 2. ip 
varies under the index shift a : n i— > n + 1. To avoid confusion, we fix the rule 
for indexing as follows: Once the points Qi, ■ ■ ■ ,Qn determined, we never 
change their induces. Alternatively, we define 

<pc(aX(y)) = <p c (X(y)) + [P-Q 1 ), 
cpc(<T 2 X(y)) = Vc {X{y)) + [P~ Qx] + [P - Q a ], 
Vc{^X{y)) = y c {X{y)) + [P- Qx] + [P - Q 2 ] + [P- Q 3 ], 
etc. ■ ■ ■ 

ycic-^Xiy)) = Vc (X{y)) - [P - Q N ], 
^c{o-~ 2 X{y)) = tp c (X(y)) - [P - Q N ] - [P - Q N -x], 
etc. ■ ■ ■ 

This particular arrangement is appropriate for our further discussion. 
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Corollary 2.13 Let 0(A(y)) be the general divisor defined by ipc(X(y)) = 
[Q(X(y)) + (N - 1) • P] (section BO). The divisor V x in pHT]) satisfies V x = 
*(*-*X(y)). 

Proof. By (j2~TTj) . we obtain 

[V x ] = [d(X(y)) -Qx Q N _x + (N - 1) • P] 

- [fKa- 1 ^)) - Qx Q N _x -Q N + N-P]. 

By the equation (y) = Q\ + - ■ - + Qn~- N ■ P, we conclude [Dx] = [X)(cr~ 1 X(y))]. 
Because T>x and ^(a^ 1 (X (y))) are general, positive and of degree g, it follows 
that Vx = d(o-- 1 X(y)). U 

As a conclusion of the corollary, we have 

(fli/fljv) - 5(A) + (A - 1)P - cKa- 1 *) - Qi Qjv^l (2.22) 

3 Tau function solution of rdKP 

Due to the linearisation theorem 12.121 and the injectivity of the eigenvector 
mapping (theorem 12. 6[) . we could say that the rndKP equation (|1.1H1.3[) is 
"essentially solved". Moreover, in some fortunate case, we can construct the 
explicit solutions by using the method of the Riemann theta functions. 

In the rest of the article, we assume that Qx = Q2, ■ ■ ■ = Qn(= Q)- Equiv- 
alently, the rndKP reduces to the rdKP equation. (See the paragraph after 
remark T2.2h 

Recall that we have assumed that g.c.d.(M + K,N) = 1 in the previous 
section. The assumption is valid also in this section. 

3.1 construction of tau functions 

We construct a theta functional solution of rdKP equation. As in the previous 
section, X t — X t (y) denotes a square matrix defined by (|2.2[) . 

Let C be the (smooth) spectral curve associated with A t . Fix a symplectic 
basis ax, ■ ■ ■ , a g ; Px, ■ ■ ■ , Pg of C and the normalised holomorphic differential 
u>x, ■ ■ ■ , u g such that J a ujj = Sij. The g x g matrix fl := (J is called 

the period matrix of C . For a fixed point po £ C, the Abel-Jacobi mapping 
A : Div(C) -> C 9 /(Z 9 + flZ 9 ) is a homomorphism defined by: 

E^-E^ - £(/>!,■■■ J> 9 )-E(J*W--,JJ 

Let us consider the universal covering 7r : il — + C and fix an inclusion l : 
C iX. For simplicity, we use the symbols 'V and "t" to express the derived 
maps Div(U) — » Div(C) and Div(C) =— > Div(il) respectively. Naturally, there 
exists a continuous lift A : Div(il) — > C 9 such that A o i(po) = 0. For the 
projection p : C 9 -> C g /{Z g + VOJ>), it follows that p o 2 = A o tt. 
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Now we should define the lifted divisors D(aX t ),D(^KXt),D(^-MX t ) € 
Div 9 (il). For fixed t € Z, assume that some lifted positive divisor D(X t ) G 
Div 9 (il) with Tr(D(X t )) — d(X t ) is specified. First, there uniquely exists a 
positive divisor D(aX t ) G Div s (il) such that: 

A(!D(aX t ))=A(p(X t ) + iP-tQ), n(J)(aX t )) = U(aX t ). (3.1) 

We will consider D(aX t ) as the appropriately lifted divisor of 5(<jX{). To choose 
appropriate D(fixlf) and S(/x_M-^t), we have to consider the compatibility: 
(hk) M + {h-m) K = id. On the Picard group on C, this reflects the equation 

[(M + iT) • P - A - A! A M -i -B q -B x = [(s)] = 0. 

f Theorem l2 . 1 21 (II) ) . Therefore, we can choose (M+iT) points kAq, • • • , kAm-i, 
kBq, • • • , kBk-i € il such that 

A((M +if)-tP-K4 «^m-i - = 0, (3.2) 

and w(KAj) — Aj, n(nBj) = Bj. We now define the two divisors D(/j,KXt), 
S(/iiwrX t ) e Div^il) by the formulas: 

= A(D{X t ) + iP- kA 3 ), n(V( PK X t )) = d(fx K X t ), (3.3) 

A{D{^ M X t )) = 2(S(X t ) + lP - nBi), ir{p(ji- M X t )) = D(^ M X t ), 

(3.4) 

where t = -(j + 1)/^ (mod M), t = -iM (mod K). 

Let r* be a holomorphic function over il defined by the formula: 

T t {jp)=e(A{Z{X t )-p-iA}), peii, (3.5) 

where #(•) = 9(»;Q) is the Riemann theta function and A £ div 9- (C) is the 
theta characteristic divisor of C ([5], Chap. II, cor. 3.11). To avoid cumbersome 
notations, we often omit the letters "A", V and use a simpler expression 
T l {p) — 9(D(X t ) — p — A), when there is no confusion possible. 

Although being defined over il, r*(p) is considered to be a multi- valued 
holomorphic function over C. By the Riemann vanishing theorem ([5], Chap. II, 
thin. 3.11), the zero divisor of r*(p) corresponds with i>(X t ). 

Let r|(p) := 9(®(aX t ) —p— A), rL(p) := 9{®(a- x X t ) -p-A). Then, by 
theorem 12. 121 the function: (a := <r _1 ) 

^ v = t'(p)-t1 + *(p) = 6(D(X t ) -p-A)- ej^KoXt) -p-A) 

satisfies [(the zeros of denominator)] = [(the zeros of numerator)] G Pic 2s (C) 
and therefore it is a single-valued and meromorphic function over C . 
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Consider an eigenvector X t (y) : = x\ ■ , (g\ = gj(x, y) = g\ (p)). 

By virtue of the equation (gt/g^) = D(X t ) + (N - 1)P - Q(aX t ) - (N - l)Q 
(|2.11[) . we derive the following equation from Liouville's theorem: 



gUp)-9i +K (p) 

Due to (|3.6[) . we can calculate some special values of ^> t (p) 



$*(p) = c x y l K ;\ yN t+K KP \ , c : constant. (3.6) 



Lemma 3.1 If g.c.d.(M + K, N) = 1, we have (i) #*(P) = c, (ii) **(Q) 

t-(M-l)K 



C X 



r *-(M-l)if ' 



Proof. Because (st^, . • • , g^ K ) = R t -(M-i)K ■ (g\, ■ ■ ■ ,9n), we have 



t . ( rt-{M-X)K t t] 



^•(/r (M - i) v 1+5 i) 



By corollary 12. 21 12.31 and proposition ^. 4[ we easily obtain the desired result. I 

Because 0(1) (X) - lQ - A) = 9(D(X) + (lP - lQ) - iP - A) = 6>(£(erX) - 
iP — A), it follows that 



**(Q) = * £ + (P), where *1 (p) 



r*(p).r^(p)- 



Then lemma O implies l\- [M ^ )K ^\(P) = ^ M ~ m ^l\P). 

Repeating the same arguments with ty+(p), we derive l\ ~^ M ^ 1 ' K ^ t ++ (P) = 
/J _(M_1)ir *^(P), and inductively we have 

Let := (P) (n "+"s). Finally we obtain the equations ** l+A r = 4^ 

and = d, where the number d does not depend on n. 

Next, consider the following single- valued meromorphic function over C: 

~t, \ = t1~ M (p)-^(p) = OWn-MSXt) - p - A) ■ 9(D(X t ) -p-A) 
[Ph T t- M (p)-Ti(p) 9(®(^ M X t )-p-A)-6(®(aX t )-p-Ay 

Using (|2.1ip and Liouville's theorem, we derive the following expression: 

9i IP) ■ 9n(p) 
which allows us to compute some special values of 
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Lemma 3.2 //g.c.d.(M + K,N) — 1, we have (i) $*(P) = d, (ii) $*(Q) = 

Trt-KM 

Proof. By (g\- M ', . . . , 5 ^ M ) = Lt~KM ■ {g{, . . . , g^), it follows that 

By virtue of corollaries 12. 21 I2TB1 and proposition ^. 41 we obtain the desired result. 
■ 

Due to $*(Q) = & + (P) and lemma [3^1 we have Vf $V(P) = V N^ t ( p )> 
which implies 

vz*\P)=v 1 -$%(p)=Vf$% + (P)=Vf$% ++ (p) = • ■ • , Vn =v?r KM ■ 

Let $^ := . + (P) (w "+"s). Therefore we obtain $* l+JV = $^ and V~& n = 
d! , where the number d! does not depend on n. 

Define := rl (tP), r* := r*(.P), r{ := r|(tP),-- ■ ,r* := r|+... + (tP) (n 
"+"s). By the above arguments, J* and V£ have following expressions: 

t+(M-l)K t+MK t+{K-l)M t+KM 

™ ~ " t+(M-l)K *+MK' ™ ~~ t+(«--l)iltf t+KM ' K J 



3.2 solution of rdKP 

For g-dimensional vectors a and fo, (a, 6) denotes a T b S C. 

Due to the periodicity OfV^Xt) = there exist integer vectors n, m S 

Z ff such that A(N(lP — iQ)) = n + ilm. Considering the definition of the 
Riemann theta function (see [5], §11.1, for example), we have 

r^ +N = x exp(— 2\^ln ■ (m, z) - \J~^\-k ■ (m, 11m)), 
where 2 = A(D(a n X t ) — iP — A). By ([3"%]) . we have 

t+(M-l)K t+MK 
Jtjt . . rt _ jJV 30 " T Af 

iiJa Jjv-a x t +(M-i)/r 

= d A 'x exp(-2\/^l tt • (m, A(lP - kAj))), (3.9) 

t+(K-l)M t+KM 

VtVj ■ ■ ■ = d' N x N 



t+{K-X)M t+KM 
' T N 



= d' N x exp(-2V^l7r ■ (m, 2(tP - kB,))), (3.10) 

where t = -jK (mod M) and * ee -iM (mod AT). Recall IL 7 ™ +M = Yln^n 
arm Ilri — Iln Ki> which imply that d depends on i (mod M) and that 

d' depends on t (mod X). Finally we obtain the conclusion: 
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Theorem 3.3 On condition that g.c.d.(M + K,N) = 1, (fXMQ3|) solves the 
rdKP equation (jj.il U.Jf wii/i constraint 



4 The general cases 

In the previous sections, we have assumed that M + K and N are relatively 
prime. This is the time to delete the assumption and discuss the general cases. 

Unfortunately, the method which we have established in this paper cannot 
be applied to the general cases. For example, when M = K = 1, N = 2, the 
defining polynomial of the spectral curve (Section O is 



where U x = + V{V£, U 2 = Vfl{ + 7 2 *i|, U 3 = I{I\VlV$. Of course, these 
Ui,U2 and U 3 are conserved quantities of the discrete reduced KP system (|l.ip . 
(|1.2p . However, there exists another hidden independent conserved quantity of 
the system. In fact, l\ + l\ + V* + V 2 l is invariant under the time evolution 
1 1— > t + 1 and is independent from U\, U 2 and U3. This means that the spectral 
curve fails to reflect faithfully the information of the system. 

Therefore, we should construct a new method for the general cases. Now we 
prove that every reduced KP equations can be traced to the case g.c.d.(M + 
K,N) = 1. Denote by KPm,k,« the reduced discrete KP equation (|1.1H1.2[) 
associated with the positive integers M, K and N. 



Proposition 4.1 (i) Suppose K < M. Define the initial values 7°, . . . , '■= 
( + o(C), (C —>■ 00, Vrt) for some complex parameter If {1^, Vfyn,t& is a 
solution ofKPM.K,N, then the sequence {X^, l^} n gz,teA converges to a solution 
ofKPM-K,K,N when ( -> 00. 

Proof. Order the elements of A as A = {ti < t 2 < ^3 < • ••}• Note that 
ts-M+i< =t s — M and 

t s -t t =0 (mod if) 1 = (mod if), (•.• t s - * s _i = 1, or K + 1). 

Especially, we have t s — t s -K = kK => t s — K, t s — 2K, . . . , t s — (k — 1)K £ A. 
To prove the statement, it is sufficient to say 



Ui = U 2 



u N . 



det (X t (y) - xE) = y 2 - y{2x + U x ) + x 2 - U 2 x + U 3 , 




Note that K < 



M A ^ and M < K & S ^ 0. 
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or equivalently, 

* = %3 + v*- kK - Vl-i + oQ.) 

yt Zr M Vt kK (1 , Q(1)) - (k&Z >0 ;t,t-kK G A). (4.1) 

By ([OHOjl and remark HJJ we have 

rt-M _ , , „,« , w , / J* = C + o(C), (Vn) 



C + o(C), (Vn 



V2 = -K*-*+o(l), (Vn) 



In our situation, it follows that t g A =*> F ? * = A ' + o(l), (Vn). Using 
(|1.1H1.2[) again, we can conclude (|4.ip soon. ■ 

Similarly, we have: 

Proposition 4.1 (ii) Suppose M < K. Define the initial values V°, . . . , V^ 1-1 
:= C + °(C)> (C oo,Vn) for some complex parameter £. if {7^, V^}„,tez is a 
solution of¥J?M,K,N, then the sequence V^}„gz i tgH converges to a solution 
of KP m,k-m,n when ( -> oo. ■ 



Example 4.2 TTie reduced discrete KP equation with M = K = 1, N = 2 can 
be traced to M = 2, K = 1, N = 2. 

W = i? = (J J), Rl = (4 1 ), and X, := 

L\RiRq. The defining function of the spectral curve is 

det (Xx(y) - S £) = -y 3 + y 2 (C 2 + E7i) - y{x(2C + E/ 4 ) + C'^i + ^3} 

where U4, = l\ + 7| + ^i. 1 + ^2 ■ Ao£e £/ia£ C/4 is i/ie hidden conserved quantity 
o/KPii : 2- if {-^n>Ki}rM * s a solution of ¥J?2,i,2, the sequence 

lim 7*, lim J 3 , lim 7^, . . . ; lim V^, lim K 5 , Hm U„ 5 , . . . 

£ — >oc £ — >oo (, — > 00 c, — >oo c, — >oo £ — >oo 

is a solution o/KPx,i,2- 
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A Proof of lemma 12.91 

Let S*, R* and L* be the matrices denned by (|2.19H2.2ip . We first calculate 
the coefficients of the equation (|2.16l42"TT7f : 

a>i,o ■ 9% + a<i,i ■ 9i+\ H h Oi,M+K ■ 9i+m+k = x ■ 9u 9i+N = V9i, 

which is equivalent to the formula X(y) ■ (g%, . ..,pj\r) = (<?i, ■ ■ • , fljv)- 
Denote the vector (gi,.. . , ffjv) by (<7i)j simply. By equation (|2.2[) : 

^t(y) : = it-(i<r-i)M(y) • • • Lt-2M(y)Lt-M{y)L t {y)x 

Rt{y)Rt-K(y)Rt-2K(y) ■ ■ ■ Rt-{M-x)K{y), 

we have 

X t ■ (ffi)i = ^t-(if-l)M • ■ • L t ^ M L t RtRt-K ■ ■ ■ Rt-{M-X)K ' (.9'0'i 

= Lt-{K-1)M ' ' ' Lt-AlLtRtRt-K ■ ■ ■ Rt-(M-2)K ' (4 <M ^^ffj + 9i+l)i 

t t d d / r t-(Af-2)K r r t-(Af-l)K . i 

= L t-{K-1)M " ••-^t-Kf ■ ••Kt-(M-3)iC " l/i i-*t 

. r,-t-(M-l)K . ~i \ 

Let A; be the set of sequences of letters s and m of length i. For example, 
Ao = 0, X\ = {s, m}, = {ss, sm, ms, mm}, ■ ■ ■ . Denote X := U;A?;. Consider 
a map (■) : X — > C defined by: (s) := 1, (m) := j* - ^ -1 ^ anc [ 

(»X):=^<X>, W := /^ M -'>* • ( X ), X e (1< I < Af), 

where er : FJ 7| • FJ i— > J| • J| Vf +1 is the index shift operator. 

Fhis notation allows us to express the vector X t • (gi)i- 
Lemma A.l Let 

Xi.k ■— {x S Xi I XTie number of 's' contained in \ is k}, 
and a,i k ■= 5~], c v (x) ■ T/iera 

' (<7i)j = ( fl i,0 ' 9% + QSi.l • + • • • + (li,M+K ' 9i+M+K)i- 

Proof. Let S; := R t -( M -l)K(y) Q = 1,2,...,M), S Af+i := L t _ {l _ 1)M (I = 
1,2,..., K) and 6 := (, = 1,2,...,M), ^ := V?"^ 1 ^ (J = 

1,2,..., K). Assume 

• • • Si • (g i ) l = (J2 k J2x£X Lk ( x ) ' 9t+k)i- 
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Then we have 

■ ■ ■ Si ■ (<7j), 
= • (Efc SxeA-i.fc M ' 5i+fe)j 

= (6+1 ' {Efc ExeAr,,* (x) ' 9i+k} + ^{Efc ExeA-i.^ (x) ' 9t+k})i 
= (Efc Exe^.fc ( m x) ' ft+fc + Efc T, xe x lik ( s x) ' ffi+fe+i)i 
= (Efe E x 'e^ 1+1 , fc (x') ■ 9i+k)i- 

By induction, we obtain X t • (fir,), = S m +a'S m +a'-i ■ • • S x • (ft)* = (a i)0 ■ ft + 

• ft+l + • • • + a^M+K ■ 9i+M+K)i- • 

In particular, we have ao/+A = (ss • ■ ■ s) = 1. 

The following lemma is given for use in the calculations further below. 
Lemma A. 2 Let x G %i,k- Then (xm) = (xs) • 7*7^ ^ 

Proof. Let 7~ := 7* ' M ^ K and assume (x'm) = (x' s ) 1 ^1/, f° r x' G Afj-i^ 
(/ - 1 > h). If x = sx' (y G <*i-i,jfe-i), we have (xm) = (sx'm) = a((x'm)) = 
<K(x'*) ■ V+fc-i) = (*x's) • = (Xs) -j; +k - If X = rn X ' (x' G *,_!,*), 
(Xm) = (mx'm) = $ • (x'm) = £; • (x's) • Ir +k = (mx's) • I i+k = ( X s) ■ I i+k . By 
induction we obtain the desired result. ■ 



A. 0.1 calculation of the determinant of 5* 

By lemma [OJ it follows that E 1 — x — oi,o = x - J2 x ex M +K,a (x)> E k+i = 
-ai,k = -E x ex M+K , k (x)- Then > 

detS* = (-1)^+1 • E X = (-l) M+K (T, x ex M+K>0 (x) - x) 

= ^M+K^t-iK-DM yt-Mytjtjt-K . . . jt—(M—l)K _ ^ 

= (-l) M+K {U 1 -x). 



A. 0.2 calculation of the determinant of R* 

Let Wk := ExeAf M x i fc (X s ) ancl ^-l ; — 0- By lemma IAT21 it follows that 

"i,k = E xe ^ M+K _ 1 , fe _ 1 (xs) + E xe * M+ K-i,* (xm) 

- Exg^m+^-x^-x (X s > + E x ^„ +K _ lti (xs) • = tu*-i + u; fc ■ 
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Let z n :— I 1 I 2 ■ ■ ■ I n . By the cofactor expansion w.r.t.the (M + K)-th row, 
we obtain: 

detiT = (-l) M+K+1 {Ei - E 2 z l + E a z 2 

+ { — J- J &M+KZM+K-1 + {— J-J ZM+ifj 

= (-1) M+A ' +1 {2: - oi j0 + 01,1*1 - ai,a«a + • • ■ 

+ (-l) M + K a ltM+ K-i ■ zm+k-i + {-1) M+K+1 z m+ k} 



= (-l) M+K+1 {x - W I 1 + W Zi - W\ZiI 2 + WiZ 2 - W2Z 2 I [ 



:i 



+ (-l) M+ ^M+K-l^M+K-l/ M+K + (-1) M+K+1 Z M+ k} 



= ( _ 1) M+X+1. X 



A. 0.3 calculation of the determinant of L* 

Recall the alternative form of the matrix X t (y) (cf. Q2.5[l ): 

X t (y) = Rt-MK{y)R t -(M+i)K{y) • ■•Rt-{2M-i)K (y)x 

Lt-(2N-i)Af(y) ■ • • L t -( K+ i) M {y)Lt-KM{y), 

in which the rightmost matrix L t -KM{y) is essential. The formula detL* = 
(— \} M + K + 1 -x comes about through arguments similar to those in section lA.0.21 
concerning (|2.5|) . ■ 
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